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ABSTRACT 

We investigate the dependence of the epoch of reionization on the properties of cosmo- 
logical perturbations, in the context of cosmologies permitted by WMAP. We compute 
the redshift of reionization using a simple model based on the Press-Schechter approxi- 
mation. For a power-law initial spectrum we estimate that reionization is likely to occur 
at a redshift z ro i on = 17^7 , consistent with the WMAP determination based on the 
temperature-polarization cross power spectrum. We estimate the delay in reionization 
if there is a negative running of the spectral index, as weakly indicated by WMAP. We 
then investigate the dependence of the reionization redshift on the nature of the initial 
perturbations. We consider chi-squared probability distribution functions with vari- 
ous degrees of freedom, motivated both by non-standard inflationary scenarios and by 
defect models. We find that in these models reionization is likely occur much earlier, 
and to be a slower process, than in the case of initial gaussian fluctuations. We also 
consider a hybrid model in which cosmic strings make an important contribution to 
the seed fluctuations on scales relevant for reionization. We find that in order for that 
model to agree with the latest WMAP results, the string contribution to the matter 
power spectrum on the standard 8/i -1 Mpc scale is likely to be at most at the level of 
one percent, which imposes tight constraints on the value of the string mass per unit 
length. 

Key words: cosmology: theory — cosmic microwave background 



1 INTRODUCTION 

The reionization history of the Universe gives important in- 
sight into the epoch of structure formation in the Universe, 
with the initial reionization believed to have been caused 
by light from the first generation of massive stars (Couch- 
man & Rees 1986; Cen & Ostriker 1992; Tegmark, Silk & 
Blanchard 1994; Fukugita & Kawasaki 1994). Over recent 
years, increasingly sophisticated techniques have been em- 
ployed to estimate the reionization epoch in cosmological 
models (Haiman & Loeb 1997; Miralda-Escude, Haehnelt & 
Rees 2000; Cen 2002; Wyithe & Loeb 2003; Somerville, Bul- 
lock & Livio 2003; Fukugita & Kawasaki 2003; see Barkana 
& Loeb 2001 for a review). 

Two important observations have begun to shed light 
on the reionization epoch, while painting a somewhat con- 
tradictory picture. The detection of Gunn-Peterson troughs 
(Gunn & Peterson 1965) in the absorption spectra of distant 
quasars suggest a very late reionization for the Universe, at 
a redshift z ~ 6 (Becker et al. 2001; Fan et al. 2003; White et 
al. 2003). By contrast, the WMAP satellite's measurement 
of temperature-polarization correlations in the microwave 



background strongly suggest a higher optical depth from 
reionized electrons, implying that reionization took place 
at a higher redshift in the range from z = 11 to z = 30 
(Kogut et al. 2003). It has been argued (Hui & Haiman 
2003) that such a high redshift of reionization would lead to 
an IGM temperature at z ~ 3 which is lower than observa- 
tions permit. In combination, these results suggest that the 
full history of ionization might be quite complicated (Hui 
& Haiman 2003; Ciardi, Ferrara & White 2003; Cen 2003) 
with more than one source of ionizing radiation. 

Until recently, one of the key uncertainties in determin- 
ing the reionization epoch was the choice of cosmological 
model, which determines the time of initial structure for- 
mation leading to the production of ionizing radiation. For 
example, Liddle & Lyth (1995) carried out the first survey 
of cosmological parameter space and found a wide range 
of predictions for different models compatible with then- 
existing observations. These uncertainties have largely been 
removed by the establishment of a Standard Cosmological 
Model based on spatially-flat low-density Universe with in- 
flationary perturbations. However, there remain significant 
uncertainties in the reionization epoch from the nature of 
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the seed fluctuations. We will study two effects, one being 
a change in the short-scale power spectrum (due either to a 
running spectral index from inflation or from a topological 
defect contribution), and the other the effect of primordial 
non-gaussianity. 

It is not our intention in this paper to provide a de- 
tailed modelling of the reionization history, but rather to 
introduce a highly-simplified model for reionization which 
allows us to assess how the epoch of reionization is modi- 
fied if the cosmological model is altered. Our model is based 
on the Press-Schechter approach, following a strategy in- 
troduced by Tegmark et al. (1994). Despite its simplic- 
ity, it makes predictions for the reionization epoch which 
are in broad agreement with those from deeper investiga- 
tions using iV-body/SPH simulations (Cen 2002, 2003; Cia- 
rdi et al. 2003) or more detailed semi-analytical modelling 
(Fukugita & Kawasaki 1994, 2003; Wyithe & Loeb 2003). 
For the purpose of surveying cosmological models, we be- 
lieve that our simple approach is well justified, and if accu- 
rate calculations become feasible for the standard cosmology 
in future, our results will enable a scaling to other cosmolo- 
gies. 



2 THE REDSHIFT OF REIONIZATION 

In the context of Cold Dark Matter Models the first stars 
form in low mass halos (~ 1O 6 M0) at high redshift. The pre- 
cise computation of the reionization redshift requires knowl- 
edge of the fraction / rc ion of baryons that must be bound 
in those halos in order to bring the fraction of ionized gas 
close to 100%. This parameter is highly uncertain since it 
depends on several astrophysical quantities, some of which 
are poorly known. Here we shall use the estimates of / rc ion 
derived by Tegmark et al. (1994). They computed the frac- 
tion of the intergalactic medium (IGM) that is ionized as a 
product of several factors, namely: the fraction of baryons 
in non-linear structures, the number of UV photons emitted 
into the IGM per proton in non-linear structures and the net 
ionizations per emitted UV photon. They concluded that the 
most reasonable estimate for the collapsed fraction required 
to induce complete reionization is / rc i on ~ 8 x 10~ 3 . This 
estimate seems to be consistent with recent more detailed 
investigations using large cosmological iV-body/SPH sim- 
ulations, where radiative transfer calculations were carried 
out by embedding massive Population III stars in gas clouds 
(Yoshida, Sokasian & Hernquist 2003). Tegmark et al. (1994) 
also computed conservative lower (/ re ion ~ 4x 10~ 5 ) and up- 
per (/roion ~ 0.8) limits for / rc ion (in their terminology re- 
ferred to as 'optimistic' and 'pessimistic') reflecting the large 
uncertainties associated with this parameter. We shall use 
these estimates in the following sections in order to estimate 
the redshift of reionization z ro ion for various models. 

Despite the large uncertainties associated with the de- 
termination of the collapsed baryon fraction necessary to 
reionize the Universe, in the context of models with initial 
gaussian fluctuations this fraction is exponentially sensitive 
to the dispersion of the density field and hence to the red- 
shift. This means that a large uncertainty in / re ion translates 
into a smaller uncertainty in the redshift of reionization, 
Zroion- However, we will see that this may not be the case 
when we consider non-gaussian fluctuations. 



3 THE MASS FRACTION 

We shall use the Press-Schechter approximation (Press & 
Schechter 1974) to compute the mass fraction, /(> M) 
associated with collapsed objects with mass larger than a 
given mass threshold M. This was originally proposed in 
the context of initial gaussian density perturbations and was 
much later generalized to accommodate non-gaussian initial 
conditions (Chiu, Ostriker & Strauss 1998). This general- 
ization of the original Press-Schechter approximation has 
successfully reproduced the results obtained from iV-body 
simulations with non-gaussian initial conditions (Robinson 
& Baker 2000). However, it has been shown by Avelino & 
Viana (2000) that this generalization docs not adequately 
solve the cloud-in-cloud problem of Press-Schechter theory 
and that deviations from the mass function predicted by this 
approach in the rare events regime are expected (see also In- 
oue & Nagashima 2002). Nevertheless, given the large uncer- 
tainties in /reion, these deviations will have a small impact 
on our final results and we shall not consider them further 
in the following analysis. 

In this context the mass fraction, /(> M), is assumed 
to be proportional to the fraction of space in which the linear 
density contrast, smoothed on the scale M, exceeds a given 
threshold S c : 

V(S)dS. (1) 

Here V(S) is the one-point probability distribution function 
(PDF) of the linear density field S and Af is a constant 
which is computed by requiring that /(> 0) = 1, thus taking 
into account the accretion of material initially present in 
underdense regions (note that Af = 2 in the case of gaussian 
initial conditions). 

Next we need to specify the filter function used to per- 
form the smoothing, and the threshold value <5 C . We shall 
use a top-hat filter 

where M — 4ivp h R 3 /3 (here p b is the background mat- 
ter density) and take S c = 1.7, motivated by the spherical 
collapse model and iV-body simulations. For spherical col- 
lapse 5 C is almost independent of the background cosmology 
(e.g. Eke, Cole & Frenk 1996). 

We consider an initial density field with a chi-squared 
one-point probability distribution function (PDF) with n 
degrees of freedom, the PDF having been shifted so that its 
mean is zero (such a PDF becomes gaussian when n — > co). 
We further assume that the shape of the PDF is indepen- 
dent of the scale and redshift under consideration, that is 
Vr{5) = P(R,5/a(R,z))/a(R,z) is always the same func- 
tion. Although this is not expected to be precisely true in 
realistic models, in particular those motivated by topologi- 
cal defects (Avelino et al. 1998a), it is reasonable to expect 
it to be a good approximation on the range of scales which 
are probed by reionization. 

In order to compute the mass fraction /(> M,z) asso- 
ciated with collapsed objects with mass larger than a given 
mass threshold M, as a function of redshift, one needs to 
compute the dispersion of the density field 
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a 2 (R,0) = k 2 \8 k \ 2 W 2 (kR)dk 



(3) 



where \Sk\ 2 is the power spectrum. For z > 1 the dispersion 
of the density field is simply given by 

a(R, 0) 



a{R ' z) ff (n° ll n£)(i + *)' 

where the suppression factor 



g(n° m ,Q° A ) = 



2 . 5^ IT1 



(4) 



(5) 



(fiS.) 4 ' 7 - n° + (i + ng,/2)(i + n« /to) : 

accounts for the dependence of the growth of density pertur- 
bations on the cosmological parameters fi^ and 0° (Carroll, 
Press & Turner 1992). This has also been verified to be a 
good approximation in the context of generic models with 
topological defects (Avelino & de Carvalho 1999). 



4 RESULTS 

Throughout, we shall adopt a cosmological model motivated 
by the WMAP results (Bennett et al. 2003). As their mild 
evidence for a running of the spectral index is driven by 
lyman-a data, whose interpretation has proven controver- 
sial (e.g. Seljak, McDonald & Makarov 2003), we shall take 
as our base cosmology throughout their preferred model as- 
suming a power- law initial spectrum (Spergel et al. 2003). 
The parameters are a matter density = 0.29, dark en- 
ergy density Q,° A = 0.71, baryon density fig = 0.047, Hubble 
parameter h = 0.72, normalization erg = 0.9, and perturba- 
tion spectral index n s = 0.99. 

To determine the amplitude of perturbations on the 
short scales relevant to reionization, for inflationary per- 
turbations we use the transfer function from Bardeen et 
al. (1986): 



\5 k \ 2 oc k Us 



ln(l + e q) 



eo<3 




-1/4 



(6) 



(T) 



(8) 



where q = k/hT, [k] = Mpc -1 , 
e= [2.34,3.89,16.1,5.46,6.71], 
and 

r = n^hexp [-n^i + v^/n^)] , 

is the shape parameter (Sugiyama 1995). We note that this 
transfer function was shown to be a good approximation 
for wavelengths small enough to come into the horizon be- 
fore matter-radiation equality if £7b = (Weinberg 2002). 
Also, although the effect of a non-zero baryon fraction is 
more complex than a simple rescaling of the shape parame- 
ter (Weinberg 2002), the Sugiyama correction turns out to 
be reasonably accurate for the values of the cosmological 
parameters we shall consider in this letter. 

4.1 The Standard Cosmological Model 

We begin by considering the simplest case, where the ini- 
tial perturbations are gaussian. Computing the perturba- 
tions and mass fraction as described, following Tegmark 
et al. (1994) and Liddle & Lyth (1995), we find a central 



value of Zroion = 17, with the range from highest to low- 
est plausible values being from 27 down to 10. This range 
is in excellent agreement with the WMAP determination of 
the reionization epoch from the temperature-polarization 
correlation function. Based on the assumption of instanta- 
neous and complete reionization, Kogut et al. (2003) find 
Zroion = 17 ± 3 (one-sigma) based on the cross-correlation 
alone, while a global fit to WMAP and other data (Spergel 
et al. 2003) gives z rc ion = 17±4. However this model appears 
in conflict with the high-redshift quasar data, and Kogut et 
al. (2003) show that allowing more complex ionization his- 
tories can widen the allowed range considerably (see also 
Haiman & Hui 2003; Ciardi et al. 2003; Cen 2003). 

The predictions above are our fiducial result, against 
which we will be comparing variations under different as- 
sumptions about the initial perturbations. If more detailed 
calculations of the reionization epoch give a different value, 
our estimates of the fractional variation about this value in 
other cosmologies should remain valid to a good approxima- 
tion. For instance, a recent detailed analysis by Fukugita 
& Kawasaki (2003) suggests a somewhat lower value of 
Zroion = 13 as the earliest possible reionization, with the 
reionization optical depth receiving a significant contribu- 
tion before reionization is complete. 



4.2 Running spectral index 

While our main focus is to examine the effects of primor- 
dial non-gaussianities, we begin with a brief exploration of 
the effect of a possible running of the spectral index, moti- 
vated by WMAP. The WMAP preferred value dn a /d\nk = 
—0.031 (Spergel et al. 2003) is much larger than anticipated 
from simple inflation models (Kosowsky & Turner 1995; 
Copeland, Grivell & Liddle 1998), and is large enough to 
have significant effect. The WMAP running model has been 
studied in detail by Somerville et al. (2003); we will consider 
a general running. 

To determine the effect of running, we define it at the 
scale k = 0.0072 Mpc" 1 at which the WMAP running cos- 
mology has the same spectral index as the power-law fit 
above. A negative running then progressively reduces short- 
scale power, delaying reionization, while a positive running 
has the opposite effect. The shift in the reionization epoch 
is shown in Figure and is well fit (to around five per cent 
within the region plotted) by 



l + 4 



1 + *3 



1 + Uy + 138y 2 



y = dn s /dink , 



(9) 



where 'nr' indicates the fiducial model without running, and 
'r' the inclusion of running. We conclude that running can 
systematically shift the expected redshift of reionization. For 
dns/dlnk = -0.031 the shift 111 1 -f~ 2reion is about 30%, 
which is consistent with the shift found by Somerville et 
al. (2003) for the WMAP running model. Since running at 
that level is permitted by the current data, we conclude that 
uncertainty in the small-scale power spectrum contributes 
a fractional uncertainty in predicting z rc ion of at least that 
level, unless one imposes a prior restricting the effect of run- 
ning. 
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Figure 1. The shift in the expected epoch of reionization induced 
by the inclusion of running. Negative running, as mildly indicated 
by WMAP, reduces the redshift of reionization. 



4.3 Non-gaussian inflation 

We now consider two possible models of non-gaussian pri- 
mordial perturbations. The first is a model with primordial 
fluctuations which may have been generated during infla- 
tion, but featuring a chi-squared probability distribution as 
may appear in the context of some non-standard inflationary 
models (see for example Salopek 1992; Linde & Mukhanov 
1997; Peebles 1999; Bartolo, Matarrese & Riotto 2001, 2002; 
Bernardeau & Uzan 2003). 

In Figure [5] we plot the evolution of the mass fraction 
/(> 1O 6 M0), as a function of redshift z, for chi-squared 
distributions with different numbers of degrees of freedom 
(n = 1, 2, 4, 8, 16, 32, 64, 128, 256, 512, oo). The intersections 
of the mass fraction curves with the horizontal lines give the 
upper limit and best estimate for the redshift of reionization. 

The lower limit is instead given by the vertical line. 
The reason is that in this case / re ion ~ 0.8 is too large to 
be reliably predicted by Press-Schechter, so we follow Lid- 
dle & Lyth (1995) in replacing this condition by the cri- 
terion that a significant fraction of the mass has to have 
collapsed for reionization to occur. We quantify this assum- 
ing a crude lower estimate for the redshift of reionization 
given by <j(10 6 Mq, z re ion) = 1 that leads to 



1 + Zt. 



= a(10 6 M e ,0)/g(n° m ,Q A ) 



which is independent of the degree of non-gaussianity of the 
initial density field. 

We see that (except in the extreme case of the lower 
limit) the predicted redshift of reionization increases as the 
number of degrees of freedom of the chi-squared distributed 
PDF decreases, as does the uncertainty on the redshift of 
reionization. We also see that the evolution of /(> 10 6 Mq) 
with redshift is not as fast for smaller n, which indicates that 
reionization may be a slower process. This happens because 
in this case /(> 1O 6 M0) is less sensitive to the dispersion 
of the density field than in the gaussian case. 

We cannot exclude any of these models, as the lower 
estimate for the redshift is independent of the degree of 
non-gaussianity. However, the best estimate increases as the 
number of degrees of freedom reduces, becoming inconsis- 




Figure 2. The evolution of the mass fraction /(> 1O 6 M ) as a. 
function of redshift z for non-gaussian inflation, as a function of 
the number of chi-squared degrees of freedom (n = 1, 2, 4, 8, 16, 
32, 64, 128, 256, 512, oo from top to bottom). The intersection 
of mass fraction curves with the horizontal lines give the upper 
limit and best estimate for the redshift of reionization. The lower 
limit is given by the vertical line. 

tent with WMAP at low n. If in the future / le ion can be ac- 
curately determined, our results will quickly indicate where 
the limit sets in. 



4.4 Inflation plus defects 

Our second non-gaussian model is one in which both string 
and inflationary perturbations are present. Although stan- 
dard topological defect models are completely excluded as 
the sole source of perturbations in the Universe, it remains 
plausible that they may play a sub-dominant role (see for 
example Durrer, Kunz & Melchiorri (2002) and refs therein). 
In particular, the popular hybrid inflation scenario relies on 
a phase transition to end inflation, which would be predicted 
to produce defects at a sub-dominant level (Copeland et 
al. 1994; Contaldi, Hindmarsh & Magueijo 1999). As we 
shall see, reionization is a particularly powerful probe of 
the effects of defects. We will focus on inflation plus cosmic 
strings, assuming initially that the two perturbation types 
have equal amplitude at the standard 8/i _1 Mpc scale. 

A good fit to the string-seeded CDM power spectrum 
is given by (Wu et al. 2002) 



|<5 fe ||oc(0.7<z) p(?) , 
where 
p(q) = 0.9 



2.7 



(10) 



(11) 



l + (2.8g)-°- 44 

Given that the inflationary and string-induced perturbations 
are uncorrelated, we can write the combined power spectrum 
as 

IMco m bi ft ed^"l^l? + (l-«)Nl , (12) 

where a is such that inflationary and string induced per- 
turbations make an equal contribution to the dispersion of 
the density field, ag on the standard 8/i _1 Mpc scale. On 
the much smaller scales relevant for reionization, the cosmic 
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Figure 3. As Figure [5] but for a hybrid model with both 
string and inflationary perturbations. The two perturbation 
types are assumed to have equal amplitude at the standard 
8h~ ^^Mpc scale. We treat the degree of non-gaussianity of the 
combined perturbation field as a variable parametrized by n = 
1, 2, 4, 8, 16, 32, 64, 128, 256, 512, oo (the number of degrees of free- 
dom of the assumed chi-squared one-point PDF). Again, the in- 
tersection of mass fraction curves with the horizontal lines give 
the upper limit and best estimate for the redshift of reionization. 
The lower limit is given by the vertical line. 



spectrum shape, with non-gaussianity further exacerbating 
this. The predicted reionization epoch is completely incom- 
patible with the observed microwave anisotropies, easily ex- 
cluding models where the defects produce half the power at 
Sh^ 1 Mpc. We note that in this case the string perturbation 
spectrum has a erf = crg/v^ ~ 0.64 which is close to the ex- 
pected value for a string mass per unit length Gfi ~ 10 -6 ) 
(Avelino et al. 1998b; Wu et al. 2002). Assuming the best 
guess for / re ion, reionization seems to occur about a factor 
of 10 earlier in redshift in the case of the hybrid model con- 
sidered above (the exact number depending on the degree 
of non-gaussianity on the scales relevant for reionization). 
Accordingly, in order to restore viability with the data the 
string contribution has to be significantly lowered, at least 
by a factor of 10. This indicates that the string contribution 
to the matter power spectrum on the standard 8/i -1 Mpc 
would have to be at most at the level of one percent, in which 
case an upper limit to the string mass per unit length would 
be Gfj, < 10~ 7 . However we note that this result is sensitive 
to the highly-uncertain value of / ro ion- Still, we can clearly 
see that reionization imposes very strong constraints on any 
contribution from defects to the seed fluctuation spectrum. 
We expect our overall results to remain valid at some level 
in the context of other defect models, as the small-scale non- 
gaussianity and excess power are generic predictions of mod- 
els of this type (Pen, Spergel & Turok 1994). 



string induced perturbations completely dominate over the 
inflationary ones. 

The cosmic string seeded density perturbations on a 
given scale 7? can be roughly divided into a nearly gaussian 
component plus a strongly-skewed non-gaussian part gener- 
ated when the string correlation length was smaller/larger 
than 7? respectively. Using simulation results for string- 
induced CDM perturbations, Avelino, Shellard & Wu (2000) 
found that the positive side of the one-point PDF was 
well approximated by a chi-squared distribution with the 
number of degrees of freedom being an increasing function 
of scale. Although these simulations did not have enough 
dynamical range to probe the scales relevant to reioniza- 
tion, they showed clearly that on length scales smaller than 
1.5(0.^ h 2 ) Mpc the perturbations seeded by cosmic strings 
have a strong non-gaussian character (Avelino et al. 1998a). 

Our results are shown in Figure |3 Here, we separate 
the effect of the non-gaussian nature of the string perturba- 
tions from that of the power spectrum. The power spectrum 
of string-induced CDM perturbations dominates on small 
scales due to the wake-like signature of cosmic strings lead- 
ing to a small-scale power spectrum close to k~ 2 . Again the 
redshift of reionization increases as the number of degrees of 
freedom of the chi-squared distributed PDF decreases, the 
same effect as that obtained for non-gaussian inflation. As 
previously discussed we expect the results obtained for small 
n to be representative of what is expected for the defect con- 
tribution (in particular cosmic strings). 

In this case we see that the much larger amplitude of 
the cosmic string contribution to the matter power spec- 
trum (relative to the inflationary one) on scales of the order 
of 1O 6 M results in a much earlier reionization (z > 65). 
Even if the perturbations were gaussian (the lowest curve) 
we have very early reionization due to the different power 



5 CONCLUSIONS 

We have shown that in a model with initial gaussian fluc- 
tuations reionization is expected to occur at a redshift con- 
sistent with the recent WMAP results, while running of the 
spectral index within the range permitted by observations 
leads to a significant uncertainty in predicting the reioniza- 
tion redshift. Although this result does not leave much room 
for non-gaussian perturbations, we have shown that even a 
small level of non-gaussianity may have interesting conse- 
quences. In the case of an initial density field with a chi- 
squared PDF, motivated both by some non-standard mod- 
els of inflation and topological defects, reionization is likely 
to be a slower process. We also show that the reionization 
history of the Universe imposes strong constraints on the en- 
ergy scale of defects given their ability to influence it through 
the induced small-scale excess power and non-gaussianity. 
We estimate that in the case of a hybrid model with cosmic 
strings the string contribution to the total power spectrum 
on the standard 8/i -1 Mpc is likely to be at most at the level 
of one percent thus requiring a low value of the string mass 
per unit length (Gfj, < 10 -7 ). 
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